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Abstract: We stude the holographic QCD model which contains a quadratic term 
— az 2 and a logarithmic term — c \og[(zjR — z)/z IR ] with an explicit infrared cut-off 
ziR in the deformed AdSs warp factor. We investigate the heavy quark potential for 
three cases, i.e, with only quadratic correction, with both quadratic and logarithmic 
O ■ corrections and with only logarithmic correction. We solve the dilaton field and 



O 



X 



dilation potential from the Einstein equation, and investigate the corresponding beta 
function in the Giirsoy -Kiritsis-Nitti (GKN) framework. Our studies show that in 
the case with only quadratic correction, a negative a or the Andreev-Zakharov model 
is favored to fit the heavy quark potential and to produce the QCD beta-function 
at 2-loop level, however, the dilaton potential is unbounded in infrared regime. One 
interesting observing for the case of positive a, or the soft-wall AdSs model is that 
the corresponding beta-function exists an infrared fixed point. In the case with only 
logarithmic correction, the heavy quark Cornell potential can be fitted very well, 
the corresponding beta-function agrees with the QCD beta-function at 2-loop level 
reasonably well, and the dilaton potential is bounded from below in infrared. At the 
end, we propose a more compact model which has only logarithmic correction in the 
deformed warp factor and has less free parameters. 
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1. Introduction 

Quantum Chromodynamics (QCD) has been accepted as the basic theory of describ- 
ing strong interaction for more than 30 years. However, it is still a challenge to solve 
QCD in non-perturbative region where gauge interaction is strong. Recently, the 
conjecture of the gravity/gauge duality [0J has revived the hope of understanding 
QCD in strongly coupled region using string theory. The AdS/CFT duality has been 
widely used to discuss the meson spectra 0, [| |[] and dense and hot quark matter 
H H 0> H> 0- The string description of realistic QCD has not been successfully for- 
mulated yet. Many efforts are invested in searching for such a realistic description by 
using the "top-down" approach, i.e. by deriving holographic QCD from string theory 
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[[TOl , [TTJ, 12], 0, as well as by using the "bottom-up" approach, i.e. by examining 
possible holographic QCD models from experimental data and lattice results. 

In the "bottom- up" approach, the most economic way might be to search for a 
deformed AdS 5 metric O, [15L O, \L7{ IT 



211, which can describe the known 



experimental data and lattice results of QCD, e.g. hadron spectra and the heavy 
quark potential. The simplest holographic QCD model is the hard-wall AdSs model 
[ p2| , which can describe the lightest meson spectra in 80 — 90% agreement with 
the experimental data. However, the hard-wall model cannot produce the Regge 
behavior for higher excitations. It is regarded that the Regge behavior is related to 
the linear confinement. It has been suggested in Ref. [13 that a negative quadratic 
dilaton term — z 2 in the action is needed to produce the right linear Regge behavior 
of p mesons or the linear confinement. 

The most direct physical quantity related to the confinement is the heavy-quark 
potential. The lattice result which is consistent with the so called Cornell potential 
[23H has the form of 

V Q Q(R) = -^ + a str R + V . (1.1) 

Where k « 0.48, a str ~ 0.183GeV 2 and V = — 0.25GeV, the first two parameters 
can be interpreted as and QCD "string" tension, respectively. 

In order to produce linear behavior of heavy flavor potential, Andreev and Za- 
kharov in Ref.[]T5f suggested a positive quadratic term modification |TJ| in the de- 
formed warp factor of the metric, which is different from the soft-wall model in 
[fT4 |. Andreev-Zakharov model has been further studied in many other articles 
p3| , [25] , |26| , |2T|| . In Ref. [p7| , the authors found that the heavy quark potential 
from the positive quadratic model is closer to the Cornell potential than that from 
the backreaction model |T7|, which contains higher order corrections. 

It is clearly seen from the Cornell potential that the Coulomb potential dominates 
in the ultraviolet (UV) region and the linear potential dominates in the infrared (IR) 
region. It motivates people to take into account the QCD running coupling effect 
into the modified metric |T9"| , pO] [TIJ. In Ref.||21||, Pirner and Galow have proposed a 
deformed metric which resembles the QCD running coupling, and the Pirner-Galow 



metric can fit the Cornell potential reasonably well. However, as shown in Ref. [28 



the corresponding dilaton potential solved from the Einstein equation is unstable, 
and the corresponding beta function does not agree with the QCD beta function. 

The motivation of this paper is to search for a deformed AdSs metric, which 
can describe the heavy quark potential as well as the QCD function and at the 
same time can have a stable dilaton potential from the gravity side. In [13|, we have 
proposed the soft-wall Dp-Dq model, which contains a quadratic correction and a 
logarithmic correction cologz. The logarithmic dilaton correction is derived from 
the top-down method, which is general for Dp — Dq system except p = 3. As pointed 
in |TJ1, that the logarithmic term Cologz itself cannot produce confinement, while 
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a logarithmic correction with an infrared cut-off in the form of Co log(zut — z) can 
have confinement at IR. Therefore, we propose a holographic QCD model, which 
contains a quadratic term — az 2 and a logarithmic term — c log[(z//j — z)/z IR ] with 
an explicit infrared cut-off Zjr in the deformed AdSs warp factor, where we assume 
a, cq can be either positive or negative. This model is found to have the same metric 
structure of the Pirner-Galow's model |Hj in UV and IR region respectively. We 
investigate the heavy quark potential in the proposed model for three cases, i.e, with 
only quadratic correction, with both quadratic and logarithmic corrections and with 
only logarithmic correction. We solve the dilaton field and dilation potential from 
the Einstein equation, and investigate the corresponding beta function in the Giirsoy 
-Kiritsis-Nitti (GKN) |Tj| framework . Our studies show that in the case with only 
quadratic correction, the produced heavy quark potential has both Coulomb part and 
linear part and fits the Cornell potential qualitatively well. In the case with only 
logarithmic correction, the heavy quark Cornell potential can be perfectly fitted and 
the corresponding beta-function agrees with the QCD beta-function reasonably well. 
At the end, we propose a more compact model which has only logarithmic correction 
in the deformed warp factor and has less free parameters. 

The paper is organized as follows. In section II, we derive the general formulae 
for heavy quark potential in the framework of AdS/CFT, and introduce the GKN 
framework to construct the gravity dual of the 5D holographic model, and calculate 
the (3 function. In section III, we construct our holographic QCD model, which 
contains a quadratic term —az 2 and a logarithmic term — Co log[(ziR—z) / Zir] with an 
explicit infrared cut-off zir in the deformed AdSs warp factor. We fit the heavy quark 
potential in this model for three cases, i.e, with only quadratic correction, with both 
quadratic and logarithmic corrections and with only logarithmic correction. We solve 
the dilaton field and dilation potential from the Einstein equation, and investigate the 
corresponding beta function. In section IV, we propose a more compact model with 
only logarithmic correction and with less parameters. The summary and discussion 
is given in section V. 



2. The formalism 

The AdSs metric in the Euclidean space takes the form of 

7-2 

ds 2 = G s flu dX"dX u = — (dt 2 + dx 2 + dz 2 ) , (2.1) 

where G s indicates the metric in the string frame, and L is the radius of AdSs. To 
search for the possible holographic QCD models, the most economic way of breaking 
conformal invariance is to add a deformed warp factor h(z) in the metric background, 



- 3 - 



and the general metric A(z) in Euclidean space has the following form: 



ds 2 = G s ltu dX lx dX v = h<KZ \ L (dt 2 + dx 2 + dz 2 ) (2.2) 
= e 2 - 4 ^) {dt 2 + dx 2 + dz 2 ) . (2.3) 

In this context, we introduce the warp factor into the pure AdSs to break the confor- 
mal symmetry to find the Q CD-like gauge theory. We will derive the general formula 
for heavy quark potential and work out the dilaton potential which will be useful 
later. 



2.1 The heavy quark potential from AdS/CFT 

To keep the paper self-contained, we follow the standard procedure [^| to derive the 



static heavy quark potential Vqq(R) under the general metric background of 
In SU(N) gauge theory, the interaction potential for infinity massive heavy quark 
antiquark is calculated from the Wilson loop 

W[C] = ^TrPexp[i I A^dx% (2.4) 

where is the gauge field, the trace is over the fundamental representation, P 
stands for path ordering. C denotes a closed loop in spacetime, which is a rectangle 
with one direction along the time direction of length T and the other space direction 
of length R. The Wilson loop describes the creation of a QQ pair with distance R 
at some time to — and the annihilation of this pair at time t = T. For T — > oo, 
the expectation value of the Wilson loop behaves as (W(C)) oc e~ TV QQ . 

According to the holographic dictionary, the expectation value of the Wilson loop 
in four dimensions should be equal to the string partition function on the modified 
AdSs space, with the string world sheet ending on the contour C at the boundary of 
AdS 5 

(W M [C]) = Z^ rmg [C}^e- s ^ , (2.5) 
where Sng is the classical world sheet Nambu-Goto action 

SnG= 2^o~ J d WDetx afc , (2.6) 

with cr s the string tension which has dimension of GeV -2 , and Xab is the induced 
worldsheet metric with a, b the indices in the (if = t, rj 1 = x) coordinates on the 
worldsheet. Under the background (|2.3j), we can obtain the equation of motion: 



2A(z) 

e 2A(zo \ (2.7) 



v / i+F) 



- 4 - 



where zq is the maximal value of z. Following the standard procedure, one can derive 
the interquark distance R function of z 

-1 e 2A{z) 1 



R(z) = 2z / du (2.8) 

Jo z A{ ] L TIma \ 2 



. e 2A(z) 
1 I .2-4(i/z) 



The heavy quark potential can be worked out from the Nambu-Goto string action: 
V QQ(z) = -^T f'dve^z— 1 (2.9) 



7l(J s JO / 1 ( e 2A(z) \ 2 

e 2A(vz) 



It is noticed that the integral in Eq. (|2.9|) in principle include some poles, which 
induces Vqq(z) — >■ oo. The infinite energy should be extracted through certain 
regularization procedure. The divergence of Vqq(z) is related to the vacuum energy 
for two static quarks. Generally speaking, the vacuum energy of two static quarks 
will be different in various background. In our latter calculations, we will use the 
regularized Vqq', for example Eqs. ( |3.3| ) and ( 4.4|) , where the vacuum energy has 
been subtracted. 

2.2 The GKN framework of the gravity dual theory 

Motivated from finding the appropriate description of heavy quark potential from 
gravity theory side, we expect to work out a general classical gravity background. 
If the deformed AdSs metric can describe QCD phenomenology, it is natural to ask 
whether it is possible to find its dual theory from gravity side. The Giirsoy -Kiritsis- 
Nitti (GKN) [19j framework offers a systematical procedure to construct the gravity 



dual theory for a 5D holographic QCD model defined in Eq. 



In this paper, we follow the notation in Ref. |28[ to introduce the GKN framework. 
According to the GKN's framework, the noncritical string background dual to the 
QCD-like gauge theories can be described by the following action in the Einstein 
frame: 

Sm-Gravity = ^| J ^^V^ (r - U^d^ - V B ((f>)) • (2.10) 

Where R is the Ricci scalar and has the dimension \R] = - — ^—r^, k% has dimension 

length 

[k 2 ] = length 3 , is the dilaton field and is dimensioness, Vb{4>) the dilaton potential 

and has the dimension of [Vr(</>)1 = — — n. The metric in the Einstein frame is 

1 yYn length 

denoted by G® which is related by the metric in the string frame G s ^ v by the following 
relation: 

G% PO = e-*+G%(X) . (2.11) 
In this subsection, the space-time metric has Minkowski signature with the sign 
convention (— , +, +, +, +). 

ds\ = e 2A{z \-dt 2 + dx 2 + dz 2 ), (2.12) 



- 5 - 



here the warp factor A(z) is in the Einstein frame, which is related to A(z) in the 
string frame by 

e 2A( Z ) = -&H?Y& = e -& e 2A(z)_ (2 .13) 

z 2 

The two independent Einstein's equations take the following form, 



3((A'(z)y + A"{z)) = -jj(0') 2 - -e^W), (2.14) 
Q(A'(z)r = 2 -{<P'f - l -e 2A ^V B {<P). (2.15) 

Adding these equations one can obtain a formal expression for the dilaton potential: 

V B (<Kz)) = -e- 2A ^(9(A'(z)) 2 + 3A"(z)). (2.16) 

Subtracting Eq. (gl4j) from Eq. (gl5|) , one can find an important relation between 
the dilaton and the metric profile: 

{<j>>Y = \({A\z)f-A"(z)). (2.17) 

It is noticed that Eq. (|2.17|) depends on the profile A(z), which is a function of the 
deformed warp factor h(z) and the dilaton field <fi{z). The resulting second order 
differential equation for <fi(z) needs two boundary conditions, which we will obtain 
from the QCD running coupling constant once the bulk coordinate z is connected 
with the energy scale E(z). 

2.3 The running coupling and the beta function 

In the GKN framework, the scalar filed or dilaton field <fi encodes the running of 
the Yang-Mills gauge theory's coupling a. For convenience, the renormalized dilaton 
field has been defined as 

a = e^. (2.18) 

The warping of the bulk space relates the bulk coordinate z to the energy scale 
E(z) via the gravational blue-shift | 28fl . By using the radial coordinate r oc 1/z, the 



blue-shift is given by the dimensionless ratio 



E r ->oo y G tt {r) 

where G u denotes the temporal component of the metric. In the limit r — > oo, the 
space-time is asymptotically flat, and G tt {r —> oo) = —1. Hence, the blue-shift reads 

Er-nn = E r yJ-G tt {r) or equivalents = E z ^J-G tt (z). (2.20) 
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A simplified expression for the energy scale in the Einstein frame has been given in 
[ffil , which has the following form: 



E r 



e -i*(*) 



a 3 



(2.21) 
(2.22) 



If one knows the value of the coupling constant a at a given energy scale E 



E r ^ooi one can find the corresponding value of z from Eq. (|2.22|). Then at given 



value of z, one can obtain <f)(z) = log(a). In order to solve Eq. ( p,17|) , two boundary 
conditions are needed. In |28| , the authors have chosen two points of running coupling 



a(E) from PDG |0 



The /3-function has the definition of 



da 

IE' 



(2.23) 



For a 5D holographic model, its function is related to the deformed warp factor 
h(z) by 



<j>'{z) 



dA A'(z) 

As we know, the QCD /3-function at 2-loop level has the following form: 



(3(a) = —boa 2 — bxa £ 



(2.24) 



(2.25) 



with b = ±(%N t 
N c = 3 and N f 



f) , and h = ^(f Nl - (f N c - ±)N f ) 0. We choose 

25 „„j l _ 77 



3 

4. In this case, 6 = |f , and 6i 



12tt 2 



The yielded beta function in |28| does not monotonically decrease with a. We 
will show in Sec. Ill that this behavior of beta function can be improved by choosing 
different boundary conditions. 



3. The holographic QCD model with quadratic and logarith- 
mic corrections 

In [|l^], we have derived the Dp — Dq model from top-down method, and found 
that for any Dp — Dq system except p = 3, there is a general logarithmic dilaton 
background field. In order to generate the Regge behavior for the light flavor mesons, 
we have proposed the soft-wall Dp-Dq model, which contains a quadratic correction 



and a logarithmic correction cologz. As pointed in |[L9|| , that the logarithmic term 
Co log z itself cannot produce confinement, while a logarithmic correction with an 
infrared cut-off in the form of Cq log(zm ~~ z ) can have confinement at IR. Therefore, 
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we extend our soft-wall Dp — Dq model to the following form with the deformed 
warp factor as 



h(z) 



exp 



az 



c ln( 



ZlR — Z, 



(3.1) 



2 ' zjr 

The coefficients a and c can be either positive or negative. An IR cut-off zm 
explicitly sets in the metric, which has the same effect as the hard-wall model M. 



When Co = 0, a > and a < corresponds to the soft-wall model [14| and Andreev 
model, respectively. 

In Ref. ||21||, in order to mimic the QCD running coupling behavior, Pirner and 
Galow proposed the deformed warp factor 



hpaiz) 



log 



(Az) 2 



(3.2) 



This metric with asymptotically conformal symmetry in the UV and infrared slavery 
in the IR region yields a good fit to the heavy QQ-potential with A = 264 MeV 
and e = A 2 / 2 = 0.48. It is worthy of mentioning that the deformed warp factor 
hpG(z) is dominated by a quadratic term az 2 in the UV regime and a logarithmic 
term — \og(zjR — z) in the IR regime, respectively. The deformed metric in Eq.( |3.1| ) 
when taking the parameter of o = 0.08, Co = 1,zjr = 2.73GeV _1 can mimic the 
Pirner-Galow deformed metric in Eq. (|3.2| ). 

Under the background (|3~T| ), the derived heavy quark potential, after subtracted 
the vacuum energy, has the form of 



V, 



QQ 



1 L 2 1 L 2 

+ 

no s z 



7T(T S z 



dv 



h{yz] 



\ 



Hz) 

h(uz) 



v- 



C Q Z 
ZlRV 



and the distance between the quark-antiquark R has the form of 



R{z) = 2z [ dvv' 
Jo 



h(z) 



1 



h{yz) 



(3.3) 



(3.4) 



h(uz) 



In the UV limit, i.e, z — > 0, the heavy quark potential has the following simplified 
expression: 



uv 0.23L 2 0.17c L 2 
v oo i K ) = — + 



+ 



a s R (T s zi R 
(0.22c + 0.2Ac 2 - 0.22a^ 2 R ) L 2 R 



a,z 



(3.5) 



IR 



- 8 - 



It is noticed that the coefficient of the Coulomb part is solely determined by the string 
tension a S} to fit the Cornell potential in UV regime, one can get a s = 0.38GeV~ 2 . It 
is also noticed that even in UV limit, both the quadratic and logarithmic corrections 
contribute to the linear potential. If Cq > and o > 0, the contribution to the linear 
potential from the quadratic term and logarithmic term compete with each other. 

3.1 With only quadratic correction 

We firstly consider the case with only quadratic correction when Co = 0. 
3.1.1 The heavy quark potential 

The heavy quark potential as functions of quark anti-quark distance R for different 
values of a = 0.1, 0.01, —0.22, — 0.4GeV 2 is shown in Fig. || (a), and the corresponding 
relation between R and z is shown in Fig.Q (b). In the numerical calculations, we 
have chosen the AdSs radius L = lGeV^ 1 , and the Coulomb part is fixed by choosing 
the string tension a s = 0.38. 




(a) (b) 

Figure 1: (a) The heavy quark potential as functions of R and (b) the distance 
R as functions of z in the case of L = lGeV" 1 , a s = 0.38GeV~ 2 , c = 0, and 
a = 0.1, 0.01, -0.22, -0.4GeV 2 . 

The behavior of R(z) is quite different for o > and o < 0. In the case of 
a > 0, the interquark diatance R firstly increases with z and reach the maximum 
R m at certain z m , then decreases with z > z m . In the case a < 0, R diverges at some 
value of z (this point is defined as z p ), the heavy quark potential also diverges at z p . 
When c = 0, it can be estimated that z p = a/— 2/ 'a. Of course, in real QCD system, 
the quark anti-quark cannot be separated to infinity. From the experimental results, 
the linear behavior breaks around R = l.lfm = 5.5GeV _1 , which is about the string 
breaking scale f3"2"|j . 
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Figure 2: The dilaton field as a function of the the bulk coordinate z in the case of 
L = lGeV -1 , a s = 0.38GeV -2 , cq = 0, and a = — 0.22GeV 2 . The boundary conditions are 
described in Eq.( |3.6D . 



When a > 0, it is found that the quark anti-quark distance R cannot reach the 
far IR regime. The larger the a is, the smaller R m can be reached. The largest R m 
for the case of a > is around 3GeV _1 , which is around 0.8fm. The slope for the 
linear potential in the middle R regime is found to be much smaller than a s t T ~ 0.183 
in the Cornell potential. 

When a < 0, it is found that the interquark distance R can go to far IR regime. 
The slope of the linear potential increases with the absolute value of |<r|. The best 
fit of the heavy quark potential gives a = — 0.22GeV 2 . With these parameters, the 
interquark distance and heavy quark potential diverges at around z p = 3.0GeV _1 . 
However, it is noticed that in this case, the slope of the linear potential is smaller 
than the experimental value. For the case of a = —0.4, the linear part is parallel 
to the Cornell potential, however, the value of Vqq is larger than the experimental 
results. Therefore, strictly speaking, the Cornell potential is not fitted very well in 
the case with only quadratic correction. 

3.1.2 The dilaton potential and the (3 function for negative a 

From the studies of heavy quark potential, we have found that a negative a is favored. 
With the parameters L = lGeV~\ a s = 0.38GeV~ 2 , c = 0, and a = -0.22GeV 2 , 
the yielded heavy quark potential is close to the Cornell potential. In this subsection, 
we solve its gravity dual and investigate the (3 function of this model. 



To solve the dilaton field from Eq. ( 2.17 ), we need to choose two boundary 



conditions. For one of the boundary conditions, we use the value of QCD running 
coupling at 3GeV as input, i.e, a(E = 3GeV) = 0.25, which can be read from 
PHI, and solve z from Eq.(|2.22|), this gives one boundary condition <ft(z = 0.87) = 
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WIHT 1 ] 0(a) 




(a) (b) 

Figure 3: (a) The dilaton potential Vb(4>) as a function of <fi and (b) The (5 function as a 
function of a in the case of L = lGeV~\ a s = 0.38GeV~ 2 , c = 0, and a = -0.22GeV 2 . 
The boundary conditions are described in Eq. fl3.6p . 

log(0.25). It is noticed that 3GeV is about the charmonium mass which is in the IR 
region. 

For another boundary condition, we can choose the same boundary condition as 
in PS| by input the running coupling at 8GeV, which gives the boundary condition 
a(8GeV) = 0.18575. However, because the produced (3 function in |28| is not a 
monotonic function, we guess this strange behavior is induced by fixing two points of 
the running coupling. Therefore, we choose to use the derivative of the dilaton field 
at z(E = 3GeV) = 0.87, i.e, <p'(z = 0.87) as another boundary condition. Because 
we don't know the value of the <ft'(z = 0.8701), we choose it as a free parameter. 

There are two types of boundary conditions we used: 

IstBC : (p(z = 0.87) = log(0.25), <p\z = 0.87) = 0.9, 

2ndBC : <f>(z = 0.87) = log(0.25), <j>(z = 0.38) = log(0.18). (3.6) 

where 4>\z = 0.87) = 0.9 is used by the best fit of the QCD (3 function, and <f)(z = 
0.38) = log(0.18) is from the input of running coupling a(8GeV) = 0.18 at UV. 

The dilaton field <p as a function of z is shown in Fig. [| It is found that 
for the two types of boundary conditions, the solution of the dilaton field <p(z) is 
monotonically increasing with z. For the 1st type of boundary condition, <fi increases 
more quickly with z than the case with 2nd type boundary condition. 

The dilaton potential Vb(4>) as a function of <p and the /3 function as a function 
of a are shown in Fig. [3] (a) and (b), respectively. It is found that for both types of 
boundary conditions, Vb{4>) decreases with 0, the dilaton potential in the IR regime 



- 11 - 




Figure 4: The dilaton field as a function of the the bulk coordinate z in the case of 
L = lGeV -1 , a s = 0.38GeV -2 , cq = 0, and a = 0.22GeV 2 . The boundary conditions are 
described in Eq.( |3.7D . 



is not bounded from below, which might indicate a unstable vacuum. For the second 
type boundary condition, i.e, the boundary condition used in PBfl , it is found that 



the produced f3 function is not a monotonic function of coupling a. This behavior 
as we have discussed, is due to the fixing running coupling constant at two points. 
For the first type of boundary condition, the produced /3 function is monotonically 
decreasing with the coupling constant a, and it agrees reasonably well with the QCD 
(5 function, which is shown by dashed line in Fig. |](b). 

3.1.3 The dilaton potential and the function for positive a 

As a reference, we solve the gravity dual and investigate the function for the case of 
positive a, which corresponds to the KKSS model or soft-wall model. The parameters 
for the model are L = lGeV~\ a s = 0.38GeV~ 2 , c = 0, and a = 0.2GeV 2 . 
The two types of boundary conditions are 

IstBC : <p(z = 0.81) = log(0.25), <j>'(z = 0.81) = 0.9, 

2ndBC : <f>{z = 0.81) = log(0.25), <j>(z = 0.38) = log(0.18), (3.7) 



which are almost the same as Eq.( 

The dilaton field as a function of z is shown in Fig. |J It is found that for the 
two types of boundary conditions, the solution of the dilaton field <j){z) monotonically 
increases to a maximum value at z m . For the first type boundary condition, z m = 
1.4GeV _1 , and for the 2nd type boundary condition, z m = l.OGeV -1 . For both 
cases, Zm is much smaller than z IR . 

The dilaton potential Vb{4>) as a function of and the (3 function as a function 
of a are shown in Fig. [5] (a) and (b), respectively. The dilaton potential Vb(4>) 
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decreases with 0, which shows an unstable potential. It is found that the f3 function 
is very interesting in the case of a positive a. The (3 function has two fixed points 
where (3 function vanishes: 1) one is a UV fixed point, where (3 = when a = 0, 
2) another is a IR fixed point, where /3 = at a moderate strong coupling constant 
a = 0.26 for 2ndBC and a = 0.32 for IstBC, respectively. If we take a lager value of 
a, the IR fixed point will appear at a smaller coupling constant, then we can have 



the Banks-Zaks fixed point [33 



3.2 With both quadratic and logarithmic corrections 
3.2.1 The heavy quark potential 

From the Pirner-Galow metric Eq.( |3.2| ), we extract the coefficient of Co = 1. Fig. |6] 
(a) shows the heavy quark potential in the case of Co = 1, the best fitted result (black 
solid line) gives a = 0.34GeV 2 and zm = 2.54GeV~ 1 . With these parameters, the 
interquark distance and the heavy quark potential diverges at z p = 1.95GeV _1 . It is 
found that the heavy quark potential is perfectly fitted in the regime R < 0.5GeV _1 
and R > 2GeV~ 1 , however, in the regime O.SGeV" 1 < R < 2GeV~\ the fitted heavy 
quark potential is a little bit higher than the experimental data. 

The result from the Pirner-Galow model is also shown in the figure by using the 
short dashed line. It is found that the Coulomb part is in good agreement with the 
Cornell potential, the linear part is parallel to the Cornell potential, however, the 
value of Vqq is a little bit higher than the data. 
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Figure 6: (a) The best fitted heavy quark potential as a function of R in the case of cq = 1 
compared with Pirner-Galow result, the UV analytical result and the Cornell potential. In 
the case of Co = 1, the other parameters are L = lGeV" 1 , o s = 0.38GeV" 2 , a = 0.34GeV 2 
and z IR = 2.54G6V" 1 . 



The green solid line is the UV analytical result from Eq. ( |3.5|) . It is found that 
this result is in good agreement with the Coulomb part of the Cornell potential in 
the region R < 2GeV~ 1 . The UV analytical is not valid any more above 2GeV _1 . 

3.2.2 The dilaton potential and the j3 function 




Figure 7: The dilaton field as a function of the the bulk coordinate z in the case of 
Cq = 1, L = lGeV" 1 , a s = 0.38GeV~ 2 , a = 0.34GeV 2 and z IR = 2.54GeV~ 1 . The two 
types of boundary conditions are described in Eq.(pT^). 
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Figure 8: (a) The dilaton potential Vb(4>) as a function of eft and (b) the (3 function as 
a function of coupling constant a in the case of co = 1, L = lGeV" 1 , a s = 0.38GeV" 2 , 
a = 0.34GeV 2 and zjr = 2.54GeV _1 . The two types of boundary conditions are described 
in Eq.(|J). 

The dilaton field as a function of z is shown in Fig. |7] for two different type of 
boundary conditions: 

IstBC : <f)(z = 1.08) = log(0.25), <p'(z = 1.08) = 2.5, 

2ndBC : <f>(z = 1.08) = log(0.25), <f>(z = 0.42) = log(0.18). (3.8) 

Where <f>(z = 1.08) = log(0.25) is from the input of running coupling a(E = 3GeV) = 
0.25 at IR, <j)'(z = 1.08) = 2.5 is by choosing the best fit of the function, and <f>(z = 
0.42) = log(0.18) is from the input of QCD running coupling a(E = 8GeV) = 0.18 
at UV. 

It is found that for these two types of boundary conditions, the solution of the 
dilaton field <fi(z) increases monotonically with z. The difference lies in that <fi is flat 
in a rather wide region of z for 2nd type of boundary condition. 

The dilaton potential Vb{4>) as a function of <fi and the /3 function as a function 
of a are shown in Fig. || (a) and (b), respectively. It is found that for the first type 
boundary condition, Vb{4>) is not bounded from below in the IR region, however, the 
/3 function monotonically decreases with the increase of a, which qualitatively agrees 
with the behavior of QCD /5 function. For the second type boundary condition, it is 
found that the dilaton potential Vb{4>) is unstable in the IR, however, the produced j3 
function is not a monotonic function of coupling a. This behavior as we have pointed 



out in Sec. 2.3, is due to the fixing of two points of running coupling constant. 
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3.3 With only logarithmic correction 
3.3.1 The heavy quark potential 

We now consider the case with only logarithmic correction when a = 0. The best 
fitted heavy quark potential as functions of quark anti-quark distance R is shown in 
Fig. |9| (a) by using the black solid line. The results are compared with that from the 
Pirner-Galow model (short dashed line) and the experimental data (the long dashed 
line) and the UV analytical result. The best fit of the heavy quark potential gives 
c = 0.272GeV 2 and Zjr = 2.1GeV _1 . With these parameters, numerical calculations 
shows that the interquark distance R becomes divergent at z p = 1.85G6V" 1 (rough 
estimates gives z p ~ Jj^|)- It is found that the heavy quark potential can be perfectly 
fitted in the whole regime of R with only logarithmic correction. 



yGeV] 



1 - 



- 1 R[GeV 



Pirner-Galow model 
Exp 

o-=0 c = 0.272 
UV analytical 



Figure 9: (a) The heavy quark potential as functions of the distance R in the case of 
L = lGeV -1 , a s = 0.38GeV" 2 , a = and c = 0.272 and z IR = 2-lGeV- 1 . 



3.3.2 The dilaton potential and the (3 function 



The dilaton field <p as a function of z is shown in Fig. [T(] for two different type of 
boundary conditions: 

IstBC : <p(z = 0.9) = log(0.25), <j/(z = 0.9) = 1.7, 

2ndBC : <f>(z = 0.9) = log(0.25), <j>(z = 0.39) = log(0.185). (3.9) 

Where <f>(z = 0.9) = log(0.25) is from the input of the a(E = 3GeV) = 0.25 at IR, 
(f)'(z = 0.9) = 1.7 is determined by choosing the best fitting of QCD (3 function, and 
(f)(z = 0.39) = log(0.18) is from a(E = 8GeV) = 0.18 at UV. 

It is found that for these two types of boundary conditions, the solution of the 
dilaton field <fi(z) increases monotonically with z. The difference lies in that <fi is flat 
in a rather wide region of z for 2nd type of boundary condition. 
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Figure 11: (a) The dilaton potential Vg as a function of (f> and (b) the /3 function as a 
function of coupling constant a in the case of L = lGeV -1 , a s = 0.38GeV~ 2 , a = and 
cq = 0.272 and zm = 2.1GeV _1 . The boundary conditions are described in Eq. (|3.9|) . 



The dilaton potential Vb{<P) as a function of <fi and the /3 function as a function 
of a are shown in Fig. [TT] (a) and (b), respectively. It is found that for the second 
type boundary condition, the dilaton potential Vb{4>) is stable which is bounded from 
below in the IR, however, the produced (3 function is not a monotonic function of 
coupling a. This behavior as we have discussed, is due to the fixing of two points 
of running coupling constant. For the first type boundary condition, Vb(4>) is also a 
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stable potential which is more deeply bounded from below in the IR region, moreover, 
the function monotonically decreases with the increase of a, and coincides with 
the QCD P function in IR regime. 

It is noticed that in the case with only logarithmic correction, the 5D dilaton 
potential Vb{4>) has the same shape as the dilaton potential in an effective 4D QCD 
model in |53J. 

4. The more compact model with only logarithmic corrections 

From studies in previous section, it is found that the model with only logarithmic 
correction in the deformed warp factor can fit the heavy quark potential perfectly, 
which is much better than the model with only quadratic correction. It might not 
be a surprise because there are four parameters used, i.e, the deformed AdS 5 radius 
L, the string tension a, the coefficient Cq and the IR cut-off Zm, while for the model 
with only quadratic correction, there are only three parameters, i.e, the deformed 
AdSs radius L, the string tension a s , and the coefficient a. 

It should be mentioned that for the case of c = 1 in Sec. |3~!2| , five parameters 
have been used to fit the heavy quark potential, and the best fitted result is better 
than the model with only quadratic correction, but not as good as the model with 
only logarithmic correction. Remind of the results in Ref. |2"7| , White found that 
the model with only quadratic correction, which has less parameters are better than 
the backreaction model, which has more parameters to produce the heavy quark 
potential. Therefore, it is not necessarily correct that one can fit the three parameters 
in the Cornell potential with enough parameters. 

Still, we hope to improve our model with only logarithmic correction. It is found 
there are two length scales in the model, i.e, the deformed AdSs radius L and the 
IR cut-off zir. We can combine these two length scales into one, and choose the 
following metric structure: 



ds 2 = G ilv dX^dX v = e 2A{z) (dt 2 + dx 2 + dz 2 ) 
= K^m , dt 2 + dS 2 + dz i 



(4.1) 



z 



with A(z) and h(z) taking the following expressions: 



A(z) = -\og(— )--£log( 



ZlR — Z 



(4.2) 



ZlR 



Z IR 




(4.3) 
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Following the same procedure, we obtain the expression of the renormalized 
heavy quark potenital VZ e J}' in the form of 



/ 



770** Z 7TO", Z 







h{yz) 



1 c z 



\ 



1 _ ^4 ( K± 

h(yz) 



2 V* Zl R V 



(4-4) 



and the interquark distance R has the form of 

R(z) = 2z [ dv 



1 e 2A(z) 



e 2A{uz) 



(4.5) 



e 2A(z) 
,2A(yz) 



The UV limit of the heavy quark potential has expression of 



Vqq(R) 



0.23z 2 R + 0.17c ^ir + (0.22c + 0.24c 2 ,) R 



a,R 



(T, 



a. 



(4.6) 



The Coulomb part can be fitted very well with the string tension a s = 1.6GeV~ , 
which is just zj R /L 2 times a s = 0.38 in Sec. |3.3.1| . The best fit of the heavy quark 
potential gives Co = 0.272GeV 2 and zjr = 2.11GeV _1 , which are the same as those 
in Sec. |3.3.1| . The results of <p(z), Vb{4>) an d /5(a) in the compact model are almost 



the same as those in Sec. |3.3.2| , therefore, we neglect the figures in this part. 

The advantage of the compact model is that with only three parameters, we can 
produce the results of heavy quark potential and QCD (3 function as good as those 
in the model with four parameters. 



5. Conclusion and discussion 

In this paper, we study a holographic QCD model which contains a quadratic term 
— az 2 and a logarithmic term — c \og[{zjR — z)/z IR ] with an explicit infrared cut-off 
zir in the deformed AdSs warp factor. We investigate the heavy quark potential, 
solve the dual gravity with dilaton field in Giirsoy -Kiritsis-Nitti (GKN) framework, 
and study the corresponding f3 function for three cases, i.e, with only quadratic cor- 
rection, with both quadratic and logarithmic corrections and with only logarithmic 
correction. Our studies show that in the case with only quadratic correction, the 
heavy quark potential can be qualitatively fitted with a negative a, and the beta- 
function agrees with the QCD beta-function reasonably well, however, the dilaton 
potential is unbounded in infrared regime. In the case with only logarithmic correc- 
tion, the heavy quark Cornell potential can be fitted very well, the corresponding 
beta-function agrees with the QCD beta-function at 2-loop level reasonably well, 
and the dilaton potential is bounded from below in infrared. We also propose a more 
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compact model which has only logarithmic correction in the deformed warp factor, 
which can describe the heavy quark potential and QCD function very well with 
only three parameters. 

Stability analysis of the dilaton potential V B 

From our numerical studies, it is shown that the dilaton potential Vb{4>) for 
the case with only quadratic correction keeps decreasing with 0, which indicates 
the potential is unstable. For the case with only logarithmic correction, the dilaton 
potential firstly decreases with then moves upward in the IR regime, which indicates 
that the dilaton potential is stable. In the following, we analyze the stability of the 
dilaton potential for the given metric Eq . (|3 . 1| ) . 

Because monotonically increases with z, we can analyze the stability of Vb{4>) 
from Vb(z) in Eq.( |2.16| ). Substitute Eq. (|2.17|) into Eq.( |2.16| ), we can have the ex- 
pression as: 

V B (z) = -e-^ (l2(A>(z)f - ^(0'(z)) 2 ) . (5.1) 



For the metric structure Eq. ( |3.1| ) in the string frame with quadratic and logarithmic 
correction, the metric in the Einstein frame takes the explicit form of 

A(z) = -~0 - A(z), A{z) = logz + -az 2 + ^log^^, (5.2) 
3 4 2 z IR 

and the dilaton potential Eq. (|5.1[ ) becomes 

V B (z) = -4et^[(0') 2 + H'A' + 3(A') 2 }. (5.3) 

Because the two square terms in the bracket are always positive, also from our 
numerical results for the physical cases, we have 0' > 0, the only chance for Vb{z) to 
change sign in the IR regime is to have a negative derivative of A. From the explicit 
expression of A', i.e, 

z 2 2{z IR - z) 

we can read that with only quadratic correction, i.e, when cq = 0, if o > 0, A' is 
always positive. When a < 0, e.g, a = —0.22 to produce the Cornell potential, A' is 
also positive in the regime of z < zir. In the case with only logarithmic correction, 
i.e, when a = 0, we can see that in the IR regime when z — > zjr, A' — > — oo when 
c > 0, therefore, Vb{z) might change sign and become positive in the IR regime. 

Positive or negative quadratic correction? To fit the heavy quark potential 
and to produce the QCD function, a negative a, i.e, the Andreev model is favored. 



However, from our previous experience in Ref.|jT3|], to produce the Regge behavior of 
p meson, a positive a is needed, i.e, the KKSS model or soft-wall model is favored. 
One possible explanation is that different model is needed to describe the physics in 
light flavor sector and heavy flavor sector, respectively. The subtlety of the quadratic 
correction in the holographic model deserves further careful studies in the future. In 



- 20 - 



our future project, we will check whether Regge behavior can be described in the 
holographic model with only logarithmic correction. 

One interesting result in the KKSS model is that the corresponding /3 function 
exists a IR fixed point. If the fixed point exists at a strong coupling regime, the KKSS 
model might be interesting to study the unitary regime of BCS-BEC crossover in cold 
atom system. If the fixed point exists at a weak coupling regime, which might be 
interesting to study the dynamical electroweak symmetry breaking physics. 
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